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ABSTRACT 

In the perturbative approach, substructures in the lens can be reduced to 
their effect on the two perturbative fields /i and A simple generic model 
of elliptical lens with a substructure situated near the critical radius is investi- 
gated in details. Analytical expressions are derived for each perturbative field, 
and basic properties are analyzed. The power spectrum of the fields is well 
approximated by a power-law, resulting in significant tails at high frequencies. 
Another feature of the perturbation by a substructure is that the ratio of the 
power spectrum at order n of the 2 fields R n is nearly 1. The ratio R n ~ 1 
is specific to substructures, for instance an higher order distortion (n > 2) 
but with auto-similar isophotes will result in R n cx -t-. Finally, the prob- 
lem of reconstructing the perturbative field is investigated. Local field model 
are implemented and fitted to maximize image similarity in the source plane. 
The non-linear optimization is greatly facilitated, since in the perturbative 
approach the circular source solution is always known. Examples of images 
distortions in the subcritical regime due to substructures are presented, and 
analyzed for different source shapes. Provided enough images and signal is 
available, the substructure field can be identified confidently. These results 
suggests that the perturbative method is an efficient tool to estimate the 
contribution of substructures to the mass distribution of lenses. 
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1 INTRODUCTION. 

Images formed in the strong gravitational regime are very sensitive to the local variations 
of the lens deflection field. In a circular potential the caustic is reduced to a single point 
and the image of a source in a singular situation is a circle. When ellipticity is introduced, 
the caustic system is a serie of connected lines with typical diamond shape aspect. This 
system presents cusps and folds singularities as predicted by the theory of singularities. 
Ellipticity is the first morphological deviation from circular symmetry, and is a standard 
ingredient of lens models. However, higher order deviations from circular symmetry are 
significant in practice, and mostly related to substructures or interaction at close range 
between galaxies, as observed in merging effects. Not all lenses are mergers, but in general 
substructures are expected in the mass distribution of lenses (see for instance Klypin et 
al. 1999, Moore et al. 1999, Ghigna et al. 2000). The presence of substructures has several 
effects, modifying the optical depth (Horesh et al. 2005, Meneghetti et al. 2007), or altering 
the image flux (Bartelmann et al 1995, Keeton et al 2003, Mao et al 2004, Maccio et al 
2006) By analyzing the image flux for very small sources, like distant quasars it should be 
possible in principle to evaluate the contribution of substructure to the lens deflection field. 
However, in practise small sources like quasars may be quite sensitive to microlensing by 
stars in the lens galaxy, which complicates the analysis. This paper will investigate the effect 
of substructures for much larger sources, typically galaxies. In particular, the modifications 
of the image morphology due to the perturbing field of the substructures will be investigated 
in details. 

2 PERTURB ATIVE DESCRIPTION OF THE EFFECTS OF 
SUBSTRUCTURE ON ARCS. 

Small substructures with mass of about a percent of the main deflector can have a major 
influence on the shape of gravitational arcs. To be effective, the perturbator must be located 
near the critical line, and in a area of image formation (|-#| < source radius). The effect of 
substructure is illustrated in Fig. ([I]) and Fig. (HJ. In Fig. ( [Q we are in the cusp caustic 
regime of an elliptical lens, while in Fig. ( [2]) the perturbative field of a sub-substructure is 
added to the elliptical lens. The field of the substructure changes dramatically the shape of 
the images, the arc is broken in 3 images, which is a situation typical of a sub-critical regime. 
In both case the perturbative approximation is over-plotted on the image of the sources ob- 
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Figure 1. Near cusp configuration for an elliptical NFW profile without substructure. The source is circular with diameter 
Rs = 0.05-Re. The black images have been obtained by direct ray-tracing. The red lines are the contours of the circular source 
obtained using the perturbative formalism. 

tained using ray-tracing. To obtain the image contours by the perturbative method, the 
derivatives of the lensing potential are estimated at r = Re (Einstein radius), providing the 
fields, fx and 4r which are directly introduced in Eq. (12) (Alard 2007). As expected the 
perturbative method gives accurate results in the elliptical case, but also in the perturbed 
elliptical case, which suggests that the perturbative method is an efficient tool to study the 
perturbations of arcs by substructures. More detailed investigations of lens with substruc- 
tures using the perturbative method are available in Peirani et al. 2008. In the perturbative 
approach all the information on the deflection potential is in two fields, h(6) and Thus 
it is sufficient to evaluate the effect of the substructure on these two fields, and this will be 
the main goal of the present paper. 

2.1 Field perturbation induced by substructures. 

For simplicity, a spherical isothermal model will be adopted for the substructure potential. 
The perturbative approach requires the estimation of the deflection field on the critical 
circle. As in Alard (2007) we re-scale the coordinate system so that the Einstein radius of 
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Figure 2. Same source position, but in addition to the elliptical deflector a substructure with mass equal to one percent of 
the elliptical lens. The perturbator is located on the X-axis at a distance r v = I.'SRe- 




the un-perturbed lens is situated at r = 1. The substructure's parameters are the following: 
its mass within the unit circle, m p , its position in polar coordinates: radius, r p = 1 + dr, and 
position angle 9 P , (see Fig. [3]). The perturbation induced by this model of substructure is 
derived in Eq. HJ The general behavior of the function's ^ and f\ is presented in Fig.'s H] 
and [5] respectively. 
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Figure 4. The variations of the field fl for different position of the substructure. In all plots the mass of the substructure is 
normalized to m p = 1 and the substructure is supposed to be on the X-axis (8 P = 0), and at the following distances from the 
center of the coordinate system: r p = 1.25 (black), r v = 1.1 (green), r p = 0.9 (blue), r v = 0.75 (red). 



2.2 Properties of the function's fi and % 

For both function's, the amplitude of the variation is quite similar for different positions 
of the perturbator, and is of the order of the mass of the perturbator. The fx function 
is asymmetrical with respect to the sign of the dr parameter (r p = 1 + dr), while for ^ 
asymmetries are weak. The other properties of the function's are related to their steep 
behavior at the origin. In particular, Fig. [5] shows that has a steep slope at the origin. 
This slope So increases as dr decrease, with approximately, S ~ 1 + -j-, which corresponds 
to a typical scale length ~ 1 — 4-. The fi functional scale like ^ near the origin, at lowest 
order fi is quadratic and the coefficient of the quadratic term is: (l + j^) 2 , thus the typical 
scale length is also: ~ 1 — -j-. 
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Figure 5. The variations of the field ^ft-, the parameters and color conventions are identical to Fig. [4] 

2.2.1 Spectral decomposition of the function's. 

The function's f\ and % are expanded in discrete Fourier series, and the power spectrum 
Pi(n) is derived from the coefficients of the Fourier expansion: 

% = J2 n °o,n cos(n6) + /3 ,„ sin(n#) 
h = Y^n a hn cos(n6>) + (3i >n sin(n0) 

(2) 

Pi{n) = al n + (3l n i = 0,1 

i. 

The power spectrum of the two function's is well approximated by a power law. The exponent 
of this power-law depends strongly on the minimum distance of the sub structure to the 
Einstein circle (see Fig. |6]). The power-law exponent in itself is variable as a function of dr, 
but interestingly the ratio of the power law components is quite constant and close to unity 
for a wide range of dr, see Fig. [TJ 



2.3 Relation to the multipole expansion of the potential 

The two fields fi and -M are related to the multipole expansion of the perturbative potential 
ip (see Eq. 3 in Alard 2007 for the definition of ijf). This relation is interesting, since some of 
the properties of multipole expansion may be exploited in the analysis of the perturbative 
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Figure 6. Power law approximation of the power spectrum of /i and . The Exponent of the power law is plotted as a 
function of the distance of the substructure to the unit circle. 
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Figure 7. This figure presents the ratio of the power spectrum's, R ■■ 
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function's. The expansion of the perturbative potential ip reads (Kochanek 1991): 

^ = -^2 (^p- cos n6 + ^4r^ sin n6 + c n (r) r n cos n6 + d n (r) r n sin n# J (3) 
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The coefficients (a n , b n ,c n , d n ) are related to the density of the lens p by the following formula: 

fln = 27m Iq 77 Jo ~ P^ u i u ) cosnv u n+1 du dv 



c « = 2^7 Jo 27 " /r=l U ) cosm; rfM rfw 

d ^ = 2^ icT X-=i v ) sinnv ul ~ n du dv 

Using Eq's (ED and (@J, and noting that: ' d ( a ™ +c ") 



(4) 



dr 



and / 
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the fields /o 



(5) 



5. 1 Local perturbation. 

Eq. [5] shows that the multipole expansion of the potential is directly related to the harmonic 
expansion of the perturbative fields. A simple and interesting case is the perturbation of the 
potential by a point mass. There are two cases, either the point mass is inside the unit circle 
(then form Eq. 0], c n = and d n = 0) or outside (o n = and b n — 0), in both cases the 
power spectrum of fi Pi(n) is equal to the power spectrum of A substructure is a local 
perturbation, and is not too far from the point mass perturbator, this analogy explains the 
ratio close to 1 which is observed in the ratio between the component of the power spectrum 
of the fields (Fig. W). 



2.3.2 Slight isophotal deformation of the density. 

Let's consider the case of a lens with small deviations from circular symmetry. In this case, 
the total density po reads: 

p (r, 6) = F (r (1 + g(9))) ~ F(r) + rF'(r) g{6) 

With: g(0) < 1 

The perturbative density p reads: 

p(r,6)=rF\r) g{9) (6) 



By introducing the former equation in Eq. (HJ) and subsequently in Eq. (EJ), the following 
equations for the fields fi and -Q are obtained: 

h = Y,n (Pn - 9») (on cos(n6>) + /3 n sin(n6>)) 



= (Pn + On) (~Pn cos(n6) + a n sin(n#)) 



W 

Pn = i; F' 



(u)u n+2 du 
q n = F \u)u 2 ~ n du 
a « = i g(v) cos(nv)dv 
P™ = h iT 9 (v) sin(nv)dv 
The ratio of the components of the power spectrum of /i, Pi(n), and, Po(n) is: 



(7) 



i2 



PM 



Pn Qn 

Pn + q n 



(8) 



For power law density profiles, F{r) oc r 7 , and R oc (-^) 2 , which is very different from 
the nearly constant ratio observed for local perturbations. 



3 PROPERTIES OF THE PERTURBED IMAGES. 

4 PROPERTIES OF THE IMAGES FORMED BY PERTURBED LENSES. 
4.1 Caustics. 

This section will investigate the caustic system of an elliptical lens perturbed by a substruc- 
ture. The lens ellipticity is defined by the parameter rj, the substructure parameters are its 
position angle 9 P , its distance dr (see Fig. [3]) and its mass m p (see Eq. [Q. For small rj, the 
elliptical potential reads: 

4>e = F (y/(l-ri) x 2 + (l + *7) V 2 ) - F{r) - | F' (r) rcos(2#) = </> (r) + i/j E (r, 9) 

Adding the contribution of the substructure (Eq. [1]), and considering that F' (1) = 1, and 
F" (1) = if the model is isothermal, the perturbative fields are: 

/ = _a cos ( 2 m + m p (i-r p cos(e-e p) ) 

2 V ; v /l-2r p cos(fl-fl p )+r2 ^ 



de V— J ' v /i_2r p cos(0-£lp)+r2 

By introducing the former analytical model in Eq. (31) from Alard (2007) we obtain analyt- 
ical equations for the caustic lines. The resulting caustic lines are presented in Fig ([8]). The 
effect of the substructure is maximal when the substructure is aligned with the potential 
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Figure 8. Caustic lines, elliptical lens, r\ = 0.1 (black), and substructure with parameters, m p = 0.03, 8 P = 0, r p = 1.3 (blue), 
and substructure, m p = 0.03, theta p = r p = 1.3 (red). The distortion of the caustic system is maximal when the the 
position of the substructure is aligned with the axis of the elliptical potential. 

axis, and increases with decreasing dr. The effect on the caustic is only weakly dependent 
on the sign of dr. 

4.2 An illustration of image anomaly due to sub-structure: sub-critical 
regime. 

The former section shows that the caustic structure is significantly modified in the vicinity of 
the substructure. An image configuration situated near the critical line for a purely elliptical 
lens is shifted to the sub-critical regime by the substructure field. An interesting point is that 
this sub-critical regime is quite different from the regime observed for the purely elliptical 
lens. The main feature of the sub-critical regime for perturbed elliptical lens is that the image 
situated at shorter distance from the substructure is much more perturbed than the others. 
Thus by comparing the structure of the images it is generally possible to detect the anomaly 
induced by the perturbator. The most obvious features of the perturbed image will be related 
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to the properties of the function's fi and -fy (see Fig.'s S]and [5j). In particular the function's 
presents large values of the slope ( or large curvature (fi) near the origin. For smaller 
images the local slope is directly related to the image size, thus the perturbation will reduce 
the size of the image. For instance let's consider the following sub-critical configurations, (i): 
a source at mid-distance from the caustic in an elliptical potential, (xs = xq, ys — 0), and 
(ii): the same system, but with the additional deflection field of a substructure with position 
angle aligned with the potential axis. In case (i) the sub-critical regime breaks the arcs in 3 
small images, the central image (8 = 0), and 2 symmetrical images. For circular sources and 
a local linear approximation of the field the size of the images is directly proportional to the 
inverse of the local derivative of -Q. Once the impact parameter is taken into account, the 
effective field in the elliptical case is: ^ = rjsm(29) — xosin(#) (see Alard 2007, Eq. 10). 
Using the former formula for % some simple calculations shows that the ratio of the size of 
the central So image to the size of the other images S\ is: 

^ = | = l + | (10) 

Eq. ( [10]) shows that for sub-critical regime, < x < 2i], thus S i > 1 consequently, the 
central image is larger that the two other symmetrical images, which is an important feature 
of the sub-critical regime for elliptical lens. With a substructure, the slope near the origin 
of ^ is perturbed by an additional field with slope at origin m p (l + j-) (see Sec. 2.2). In 
the hypothesis of a local perturbation the 2 other images remain unperturbed. Assuming 
dr the size of the central image is now: 

So 1 oc ^ + (2rj - x ) (11) 

The scale of x in the caustic system is rj, thus Eq.'s ( [TT]) indicates that the image ratio is 
modified by: 

Sq! = Sqi X — r ^-m£ (12) 
" + dr 

Thus the perturbation is significant if at least, m p ~ rj dr. Since the usual scale of both r\ 
and dr is of about one tenth of Re, m p must be of only the order of a percent to alter very 
significantly the size of the central image. 



4.3 Analysis of images by Reconstruction of the perturbative field. 

This section will show that the image anomalies due to substructures can be reconstructed 
independantly of the source shape. The reconstruction will be illustrated by the configuration 
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presented in Sec. 4.1. In this configuration, the distortion is very obvious, because most of 
the effect is on the central image, and that the size of this image is reduced by the field 
of the substructure. However this results holds for source with circular contours only. For 
other sources the size of the images may be modified slightly, and for better accuracy a 
more general approach is required. Here we have to tackle the general problem of lensing 
potential and source reconstruction, taken together these problem are difficult and may be 
quite degenerate. The problem in itself is much simplified if we make the hypothesis that the 
field is smooth at the scale of the image and can be represented by a lower order polynomial in 
6. The knowledge of the field at the images positions allows to transfer the image contours to 
the source plane, where the different images of the sources must be identical. The constraint 
that the image must be identical in the source plane will be used to evaluate the field at the 
image position. Another constraint on lens reconstruction is that no images are produced 
in void areas (Diego et al ). This constraint is simple to implement in the perturbative 
approach, it is sufficient that \^\ > Rc, where Rc is the radius of a circular envelope to 
the source contour. 

4-3.1 Local field models. 

There are basically two kind of image models to consider, smaller images, with nearly linear 
field at the scale of the image, and the longer images produced in near caustics configurations 
. For smaller images, the model to adopt for the f field is very obvious, ^ ~ a (6 — (3 ), 
with cto a local constant to evaluate, and (3q a position angle which in practice should be 
close to the image center. Near caustics the model should be a polynomial of higher order. 
For the f\ field the situation is similar to the fo field for small images, while in near caustics 
situations the local modelisation of the fi field will require a full polynomial expansion to 
higher order (2 or 3). 

4-3.2 Fitting local field models. 

Simulated data are obtained by ray-tracing the lens model in the sub-critical configuration 
described in Sec. 4.2. The ellipticity of the lens is rj = 0.1, and the substructure parameters 
are:, rap = 0.03, rp = 1.3, dp = 0. The source position is xs = I.577, ys = 0. The fitting 
procedure is greatly simplified by the fact that in the perturbative approach the solution 
for a circular source is known (Alard 2007). For circular source the f\ field corresponds 
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to the mean radial position of the contour at a given 9, and the field is the contour 
width. Once these fields are extracted from the images data, they make a good starting 
guess for the real solution. Basically the guess will be estimated by fitting local polynomial 
expansions of 9 to the circular solution. Starting from this guess a few steps of non-linear 
fitting optimization will lead to the correct solution. The quantity to minimize to achieve an 
estimate of the solution will be a measurement of the image similarity when a re-mapping 
to the source plane is performed. A general method to estimate the image similarity is to 
compare their moments. Suppose that there are Nj images, and that their moments to order 
Ns in the source plane are identical, this provide N s( N s+3)(N!-i) constraints. Local linear 
image models have 2iVj parameters, thus even in the case of 2 images only the problem is 
already over-constrained for Ns = 2. In the case presented in Fig's ( E]), ( [TO]) and ( ITTT) . 
there are 4 images, thus moments up to the second order are sufficient to fit the data. Ray 
tracing of the sub-critical configuration presented in Section 4.2 is performed for different 
sources shape, resulting in 3 different sets of 4 images. Local linear field models are fitted to 
each of these sets. The circular solution is used as a first guess and a quantity that measures 
the distance between the image moments in the source plane is minimized using the simplex 
method. In each case both the field parameters and source shape can be recovered (See table 
[1]). The great advantage of the local field method, is that even for the higher order fields 
generated by a substructure, the local models are simple and can be recovered. Note also that 
for instance the slope of -Q gives 4 constraints, which coupled with the image positions gives 
8 constraints. In the case of an elliptical lens the Fourier expansion is of order 2 (provided 
ellipticity is not too large), which corresponds to 4 parameters only. Thus the elliptical case 
is over-constrained in the present situation. With only 2 images the elliptical lens is already 
constrained. Taking any 2 images unperturbed by the substructure, and making an elliptical 
model, it will always be impossible for this model to meet the constraint related to the 
perturbed image. 

4-3.3 Accuracy of measurements using the perturbative method. 

In Table 1 the perturbation on the local slope of the ^ function due to the substructure is 
about 7 times the mean scatter of the measurements recovered by fitting local models for the 
different sources. Translated in accuracy on the measurement of the substructure mass gives 
about 0.4 % of the main halo mass. However the scatter in the measurements comes from the 
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Table 1. Local model fitting results for different sources, first circular source with diameter Rg = 0.05-Rg> elliptical source 
with r]g = 0.5, (f — r]g)x s + (1 — T]s)y s = Rg, and Rg = 0.05Re, two circular sources with diameter Rg = 0.05-Rg and 
respective centers: (— ^-,0), (-^-,0). For comparison, in the last column of the first part of the table, the theoretical results 
of the local field values are given when no substructure is present (see the large difference on the slope of image 2). It is not 
presented for the second part, since due to symmetry of the substructure position it has no effect on the /i field. 
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Figure 9. Image obtained by ray-tracing for a circular source. The lens is elliptical, r)g = 0.1, and is perturbed by a substructure 
with parameters: m p = 0.03, theta p = 0, r p = 1.3. The source has diameter Rg = 0.05i?E, and impact parameter, (Xg = 
0.15i? B ,0). 
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Figure 10. Same as figure [9] but for an elliptical source with ellipticity r]s = 0.5 

perturbative approximation and also from the simplicity of the local (linear) modeling. Thus 
this scatter is an over-estimate of the error made in the perturbative approximation. The 
perturbative method may introduce some limitation in accuracy, but in practice, resolution 
effects and noise should be much stronger limiting factors. And more importantly, the limit 
in accuracy by the perturbative method may be a concern only for absolute measurement 
of the displacement field of the lens, if we are interested in the differential effect of the 
substructure field, then the perturbative method will be very accurate, first because the 
method is linear, and thus allows differential measurements, and second because also due to 
linearity, the field of the perturbator can be reconstructed with an accuracy that scales likes 
its mass. Thus, for the evaluation of the differential effects due to very small substructure 
the perturbative method should give accurate results, in this case, the main problem will 
be related to the un-biased statistical estimation of the perturbative fields for distributions 
without substructure. 
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Figure 11. Same as figure [9] but for two circular sources with diameter Rg = 0.05 and respective centers: (- 



,0), (^,0). 



4.4 Approximate source invariant quantities 

In general the image features are dependent upon the source shape, and the re-construction 
of the lens fields require the non-linear procedure presented in the former section. However, 
for smaller images and weakly elliptical sources, there is an approximate invariant. This 
conserved quantity is useful for nearly round sources, or for improving the first guess in 
the non linear fitting procedure. Let's consider small images with total size Qj <C 1 and an 
elliptical source, the width of the image is given by Eq. (15) in Alard (2007): 



(13) 



w = ^ s ~ {1 ~f s) ^ 

S=l- Vs cos (2 (0-0„)) 
Where rjs is the source ellipticity and 9q is the angle of the source main axis. Since the 
image is supposed to be small, we are operating in a small range of 9, and the field ^ may 
be linearized locally. Near the center of the image, ^ ~ 0, thus by taking the origin of 9 
at the image center we have: ^ ~ k9. We will also assume that the ellipticity is a small 
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number, so that by change of variable r/s = er/s } and 6 = e6. Using these new variables, 
it is possible to expand Eq. ( [TBI) in series of e, which simplifies the calculation of many 
quantities. . In particular, the image size along the orthoradial direction is obtained by the 
condition W = 0. Solving to the lowest order in e, we obtain a second order equation in 9, 
and the difference of the 2 roots gives the image size. To the lowest order in e, the image 
size in the orthoradial direction is Wt'- 



The size of the image in the other direction is approximately the size of the image in the 
radial direction near the center of the image, from Eq. ( /refwidth) to the lowest order in e, 
the radial size Wr is: 



To first order in rjs, the product S = Wt x Wr which is closely related to the image surface 
does not depend on the ellipticity of the source. This result means that in practice for small 
ellipticity (rjs <S 1), S is a constant independent of the source ellipticity. 
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